Hole spin dephasing in p-type semiconductor quantum wells 
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Hole spin dephasing time due to the D'yakonov-Perel' mechanism in p-type GaAs (100) quantum 
wells with well separated light-hole and heavy-hole bands is studied by constructing and numerically 
solving the kinetic spin Bloch equations. We include all the spin-conserving scattering such as the 
hole-phonon and the hole-nonmagnetic impurity as well as the hole-hole Coulomb scattering in our 
calculation. Different effects such as the temperature, the hole density, the impurity density and 
the Rashba coefficient on the spin dephasing are investigated in detail. We also show that the 
Coulomb scattering makes marked contribution to the spin dephasing. The spin dephasing time can 
either increase or decrease with temperature, hole/impurity density or the inclusion of the Coulomb 
scattering depending on the relative importance of the spin-orbit coupling and the scattering. It is 
also shown that due to the different spin-orbit coupling strengthes, many spin dephasing properties 
of holes are quite different from those of electrons. 

PACS numbers: 72.25.Rb, 71.10.-w, 67.57.Lm, 78.90.+t 



I. INTRODUCTION 

Much attention has been devoted to the spin degree 
of freedom of carriers in Zinc-blende semiconductors re- 
cently due to the possible application to the spintronic 
devices. 1,2,3,4 Understanding spin dephasing/relaxation 
of carriers in semiconductor quantum wells (QWs) is 
one of the most important prerequisites for the realiza- 
tion of such devices. There are many studies of spin 
dephasing/relaxation of electrons in n-type QWs where 
the spin dephasing is determined by the D'yakonov and 
Perel' (DP) mechanism. 5 Most studies are within the 
framework of single-particle approximation 6,7,8,9,10,11,12 
and the spin dephasing time (SDT) can be written as 6 
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Here T p (k) is the momentum relaxation time. Jl(k) is the 
DP term which is composed of the Dresselhaus term 13 
due to the bulk inversion asymmetry (BIA) and the 
Rashba term 14 due to the structure inversion asymme- 
try (SIA). fi(k) 2 denotes the average of Jl 2 (k) over all 
directions of k. In GaAs QWs, the Dresselhaus term is 
the leading term and fi(k) has the form: 



n x (k) = 7 k x (k 2 y - <fc 2 )) 
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in which (fc 2 ) represents the average of the operator 
— (d/dz) 2 over the electronic state of the lowest subband 
and is therefore 7r 2 /a 2 under the infinite- well-depth as- 
sumption. 7 is the Dresselhaus spin-orbit parameter. 6,15 
It is noted that Eq. (1) is valid only when |f2|T p <C 1 
and the scattering is elastic. In this limiting case, the 
angular rotation of electron spin over time t p is small 
and spin relaxation occurs as a result of a number of 



accidental small rotations. 6 This approach captures the 
lowest (first) order of the anisotropy due to the fact that 

ft(-k) = -n(k). 

It is shown recently by Wu et al. from a full many- 
body microscopic a pp roac h 16 ' 17 ' 18 ' 19 ' 20 ' 21 ' 22 ' 23 that the 
single-particle approach is inadequate in accounting for 
the spin dephasing/relaxation. The momentum depen- 
dence of the effective magnetic field (the DP term), and 
even the momentum dependence of the spin diffusion 
rate along the spacial gradient 22 or the random spin-orbit 
coupling, 24 serve as inhomogeneous broadenings. 17,18 In 
the presence of the inhomogeneous broadening, any scat- 
tering, including the carrier-carrier Coulomb scattering 
(beyond the Hartree-Fock self-energy from the Coulomb 
interaction), can cause irreversible dephasing. This 
many-body approach takes account of the inhomoge- 
neous broadening not only from different directions of 
n(k) (not only -|0(k)| and +|fi(k)|), but also from the 
modulus of the DP effective field, i.e., |fi(k)|. 25 More- 
over, this approach also takes full account of the counter 
effect of the scattering to the inhomogeneous broaden- 
ing instead of only the lowest-order elastic scattering: 
The scattering tends to drive carriers to more homoge- 
neous states and therefore suppresses the inhomogeneous 
broadening induced by the DP term. Finally, this ap- 
proach is valid even when |fl(k)|r p gl 1 and is applicable 
to systems far away from equilibrium (eg., systems with 
high spin polarization 19 and/or in the presence of high 
electric field parallel to QWs). 20,21 Using this method, 
Weng and Wu performed a systemic studies of spin de- 
phasing in n-type GaAs (100) QWs at high temperatures 
and showed that the effects beyond the single-particle 
approach Eq. (1) are dominant even for systems near 
equilibrium. 19 These effects include the many-body ef- 
fects, the inhomogeneous broadening induced spin de- 
phasing and the counter effect of the scattering to the 
inhomogeneous broadening. For small well width, the 
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calculated electron SDTs using this microscopic many- 
body theory increase with temperature and are in agree- 
ment with the experiment both qualitatively and quan- 
titatively, while the SDTs of earlier simplified treatment 
drop dramatically with temperature and are one order of 
magnitude larger than the experiment data. 19 For larger 
well width, the SDT may first increase then decrease with 
temperature. 21 These properties come from the compet- 
ing effects between the DP term and the scattering. 

Although there are extensive investigations on the spin 
relaxation/dephasing of electrons, investigations on the 
spin relaxation/dephasing of holes in p-type semiconduc- 
tor QWs are relatively limited. 26 ' 27 Nevertheless, knowl- 
edge of the spin relaxation/dephasing of holes in p-type 
QWs is very important to the assessment of the fea- 
sibility of hole-based spintronic devices. This is be- 
cause a possible way to achieve high electronic spin in- 
jection without the conductance mismatch 32 is to use 
magnetic semiconductors as spin source and most mag- 
netic semiconductors are p-type at high temperature. 33 
Very recently there are some reports on the hole spin 
relaxation/dephasing. 28 ' 29,30 ' 31 All the theoretical calcu- 
lations in these works are within the framework of the 
single-particle approximation Eq. (I). 30 ' 31 

It has been shown in electron systems that Eq. (1) is 
inadequate in accounting for the spin dephasing. More- 
over, the electronic states and spin-orbit coupling of holes 
are very different from those of electrons. 34,35 ' 36 In bulk 
material, the T-point degeneracy of the heave hole (HH) 
and the light hole (LH) makes the hole spin relaxation in 
the same order of the momentum relaxation (100 fs). 28 
This degeneracy is lifted in QWs. Under the parabolic 
approximation, the HH and LH bands can be treated in- 
dependently for QWs of small well width. Unlike the con- 
duction band where the DP term mainly comes from the 
BIA contribution in GaAs QWs, in p-type GaAs QWs, 
the SIA contribution is usually the dominant one. It is 
noted that in hole system the relation \fl\r p <C 1 is usu- 
ally unsatisfied due to the strong spin-orbit coupling and 
consequently the validity of Eq. (1) is even more ques- 
tionable. Therefore, in this paper we investigate the hole 
spin dephasing using our full many-body microscopic ap- 
proach. We calculate the SDT of the HH and LH by nu- 
merically solving the many-body spin kinetic Bloch equa- 
tions with all the scattering explicitly included. Then 
we discuss how the temperature, the hole density, the 
Coulomb scattering, the Rashba coefficient and the im- 
purity density affect the SDT. We show that the eariler 
treatment based on the single particle approximation is 
not valid in hole systems and unlike the case of electrons 
where the scattering "always" raises the SDT at low-spin 
polarization, the scattering can either enhance or sup- 
press the SDT of holes based on the relative importance 
of the Rashba term and the scattering. 

The paper is organized as follows: In Sec. II we set 
up our model and kinetic equations. Then in Sec. Ill 
we present our numerical results. We first show the time 
evolution of the spin signal in Sec. Ill A. In Sec. Ill B 



we investigate how the temperature affects the spin de- 
phasing. The Coulomb scattering, the impurity density 
and hole density dependence of the SDT are discussed 
separately in Sec. Ill B, C and D. We conclude in Sec. 
IV. In Appendix A we show the effect of the scattering 
to spin dephasing when it is much weaker than the spin- 
orbit coupling strength. In Appendix B we present a 
simplified analytical analysis of the SDT and show the 
different effects of the scattering to the spin dephasing at 
strong/weak scattering regime. 



II. KINETIC SPIN BLOCH EQUATIONS 

We start our investigation from ap-doped (100) GaAs 
QW of well width a. The growth direction is assumed to 
be along the z axis. A moderate magnetic field B is ap- 
plied along the x axis (in the Voigt configuration). Here 
we assume only the lowest subband is populated. It is 
noted that for two-dimensional hole system, the lowest 
subband is HH-likcd. By applying a suitable strain, the 
lowest subband can be LH-liked. We assume the con- 
finement is large enough so that the HH and LH bands 
are well separated and we may consider the HHs and LHs 
separately. With the DP term included, the Hamiltonian 
of the holes can be written as: 

Hx = ^ 5 ™' + + n\k)} ■ ^l} c t kff£W 

kcrcr' 

+ Hj . (3) 

Here A = LH, HH denotes the LH or HH state, a = +, — 
stands for the spin, e^x = k 2 /2m* x is the energy of hole 
with wave vector k and effective mass m* x . a are the 
Pauli matrices. The DP term is mainly from the Rashba 
term. For (100) GaAs QWs, we have 
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for HHs and 

^(k) 



+"fi!i l ky(ky — 3/c 2 )] , 

n^(k) = -2E z [^ 6l (k 2 z )k x +j! l 3 6l kp x 

+ltTk x {kl-3kl)] , 

n^ H (k) = o 
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for LHs. 35 It is seen from these equations that the mag- 
nitude of the Rashba term can be tuned by means 
of an external gate voltage which changes the elec- 
tric field E z in the sample. 38 ' 39 ' 40 - 41 7 ™ 7/ \ j 7 5 f h 
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energy gaps between the HH and the LH bands, the HH 
and the split-off bands, the LH and the split-off bands 
respectively: 
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with Ao representing the energy gap of the split-off band 
(from the T-point of the valence band). 71, 72 and 73 
are the Luttinger parameters. From Eqs. (4-9) one can 
see that the HH Rashba terms include only the cubic 
terms whereas the LH ones include both the cubic and 
the linear terms. The ratio of the cubic and the linear 
terms depends on the well width a: (k 2 ) in the linear 
terms decreases with a 2 . Futhcrmore, one can see from 
Eqs. (10) and (12) that A^; decreases faster with a than 
A/ s , which makes 7|| 6/ and -ff 1 ® 1 increase faster with a 
than 7|2 6 '- Therefore, the cubic terms weighted by 7|| 6/ 



and 7|4 6 ' increase faster with well width than the linear 
terms weighted by 752 s ' • In brief, when a is small, both 
the linear and the cubic terms are important; When a 
gets larger, the cubic terms are the leading terms. 

The interaction Hamiltonian Hi in Eq. (3) is composed 
of the hole-hole Coulomb interaction and hole-phonon 
scatering, as well as hole-impurity scattering. Their ex- 
pressions can be found in textbooks. 42,43 

We construct the many-body kinetic spin Bloch equa- 
tions by the non-equilibrium Green function method as 
follows: 



PkX.aa' PkX.aa' \coh ~r pkX.aa' \ scatt 



(13) 



with pkA.crcr' representing the single-particle density ma- 
trix elements. The diagonal elements p k A,<T<r = /kA,<x de- 
scribe the hole distribution functions of wavevector k , 

state A and spin a. The off-diagonal elements Pk.A.n = 

p k A h = pkA describe the inter-spin-band correlations 

for the spin coherence. p\x. a a'\coh describe the coher- 
ent spin precessions around the applied magnetic field 
B in the Voigt configuration, the effective magnetic field 
f2 A (k) as well as the effective magnetic field from the 
hole-hole Coulomb interaction in the Hartree-Fock ap- 
proximation and can be written as: 
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p\\,aa' \scatt in Eq. (13) denote the hole-hole Coulomb, 
hole-phonon and hole-impurity scattering. The expres- 
sions of these scattering terms and the details of solv- 
ing these many-body kinetic spin Bloch equations can be 
found in Rcf. [20]. ' 



III. NUMERICAL RESULTS 

We numerically solve the kinetic spin Bloch equations 
and obtain temporal evolution of the hole distribution 
fkx,a(t) and the spin coherence pk\(t). We include the 
hole-phonon and the hole-hole scattering throughout our 
computation. As we concentrate on the relatively high- 
temperature regime (T > 120 K), we only include the 
hole-LO-phonon scattering. The hole-impurity scatter- 
ing is included when stated. As discussed in the previous 
papers, 16 ' 17,42,44 the irreversible spin dephasing can be 
well defined by the slope of the envelope of the incoher- 



ently summed spin coherence 

Px = 52\Pkx(t)\. (16) 

k 

The material parameters of GaAs in our calculation are 
tabulated in Table I where fi LO represents the LO phonon 
frequency and Kqo (kq) is the optical (static) dielectric 
constant. 35,45 Our main results are plotted in Figs. 1- 
6. In these calculations the width of the QW is chosen 
to be 5 nm unless otherwise specified; the initial spin 
polarization P x = (N x .+ - N X ^)/(N X . + + A A: _) is 2.5 % 
with 

^A,a = ^/kA,a, (17) 
k 

representing the hole density of cr-spin band; the mag- 
netic field B = 4 T, and the Rashba coefficient 
"f™ h E z m is taken to 0.5 nm when a = 5 nm. 
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FIG. 1: Typical hole densities of spin-up and spin-down bands 
and the incoherently summed spin coherence phh vs. time t 
shown for the case of the HHs. Note that the scale of the spin 
coherence is on the right side of the figure. The dotted line 
represents the slope of the envelope of phh- 



A. Temporal evolution of the spin signal 

We first study the temporal evolution of the spin signal 
in a GaAs QW at T = 300 K. In Fig. 1 we show the 
typical evolution of the HH densities in the spin-up and - 
down bands together with the incoherently summed spin 
coherence for the total HH density Nhh,+ + Nhh.- = 
4 x 10 11 cm~ 2 and impurity density Ni = 0. It is seen 
from the figure that the hole densities in the spin-up and - 
down bands and the incoherently summed spin coherence 
oscillate due to the presence of the magnetic field. From 
the slope of the envelope of the incoherently summed spin 
coherence, one is able to deduce the SDT. 



TABLE I: Parameters used in the calculation. 
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B. Temperature dependence of the SDT 

We now turn to study the temperature dependence of 
the SDT at different impurity densities Ni. We plot the 
SDTs of the LH and the HH in Fig. 2(a) and (b) as func- 
tions of temperature. The total LH and HH densities 
Nlh and Nhh arc taken to be Nh = 4 x 10 11 cm~ 2 . 
One finds from Fig. 2(a) that for the LH, when there 
are no impurity Ni = or low impurities Ni = O.lNh, 
the SDT first decreases then increases with temperature. 
The minimum occurs at smaller temperature for higher 
impurity densities: 140 K when N = 0.1N h and 200 K 
when Ni = 0. When the impurity density Ni — Nh, the 



SDT increases with temperature monotonically. These 
temperature dependences are quite different from those 
of electrons in QWs with the same electron density and 
initial spin polarization where the SDT increases mono- 
tonically with temperature. 19 

It is noted that the property of spin dephasing is 
quite different when |fi|r p ^ 1 and \SI\t p <C 1. When 
\n\T p ^ 1, the scattering is weak in comparison to the 
DP effective field (inhomogcncous broadening) and the 
counter effect of the scattering to the inhomogeneous 
broadening is unimportant or can be ignored. In the pres- 
ence of inhomogeneous broadening, adding a new scat- 
tering provides an additional dephasing channel. 17,42,46 
This effect has been revealed in detail in Appendix A. 
Therefore, the scattering in this regime provides a spin 
dephasing channel and the increase of the temperature 
leads to a stronger scattering and consequently a faster 
spin dephasing. Moreover, the increase of the tempera- 
ture drives holes to a higher k-state, and holes experience 
a larger |fi(k)|, i.e., a stronger inhomogeneous broad- 
ening. This tends to reduce the SDT too. Therefore, 
the SDT decreases with temperature when \fl\r p 1. 
When \ fl\Tp <C 1, the scattering is strong in comparison 
to the DP term. Hence the counter effect of the scattering 
to the inhomogeneous broadening cannot be ignored any 
more. As the scattering tends to drive carriers to more 
homogeneous states in k-space, it tends to increase the 
SDT. Therefore, whether the SDT increases or decreases 
with temperature depends on the competition between 
the scattering and the DP term. It will be shown later 
that when the linear part in the DP term is dominant, the 
increase of the inhomogeneous broadening with temper- 
ature is relatively slower than that of the scattering and 
the SDT increases with temperature. Nonetheless, when 
the cubic part in the DP term is dominant, the increase 
of the inhomogcncous broadening with temperature turns 
out to be faster than the increase of the scattering and 
the SDT decreases with temperature. 

For electrons in GaAs QWs, the spin-orbit coupling 
is not very strong. \SI\t p is usually much smaller than 
1 (typically \S1\t p = 0.016 at T = 100 K, a = 15 nm, 
and the total electron density N e ^ + + N Si _ = 4 x 10 11 
cm -2 ). Therefore when the linear (cubic) term in Eq. (2) 
is dominant, the SDT of electrons increases (decreases) 
with temperature. 19,21 



TABLE II: Rashba coefficients [unit: nm/(B 2 mo)]. 





6161 


„ 6161 


„ 6161 


„ 7h7h 


_ 7h7h 




752 


753 


754 


753 


754 


a — 5 nm 


-0.193 


0.650 


0.089 


-0.080 


0.500 


a — 7 nm 


-0.156 


2.21 


0.341 


-0.330 


2.07 



Situations are more complicated for hole system due 
to the strong spin-orbit coupling. The Rashba coeffi- 
cients of the coupling in Eqs. (4-9) are listed in Table II. 

For LHs, when a = 5 nm, E z mo¥rix changes from 

—3.36 nm to —1.46 nm when the temperature changes 
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FIG. 2: SDT r of LHs [Fig. 2(a)] and HHs [Fig. 2(b)] vs. 
temperature T at different impurity densities. ■: N = 0; 
A: Ni = 0.1N h ; and •: iV; = 1.0N h . The dashed curves in 
Fig. 2(a) are the corresponding curves but at bigger well width 
a — 7 nm. The SDTs calculated from the simplified treatment 
(solid curves) and our many-body treatment (dashed curves) 
are plotted in the insets for comparison. 

from 100 K to 300 K. Here (A; 2 ) represents the average of 
fcjj. It is seen from Table II that both jf® 1 and 754 s ' are 

smaller than — Ts2 6Z "^^y • Therefore the linear terms in 

Eqs. (7) and (8) are dominant. From these coefficients, 
one may further find that when there is no impurity, the 
spin-orbit coupling for LHs is one or two orders of magni- 
tude larger than that of electrons. Consequently neither 
\fl LH \r p > 1 nor \Q, lh \t p < 1 is satisfied here. The 
value of IJl^Irp is usually slightly smaller than 1. In 
this regime, both the competing effects of the scatter- 
ing addressed above can not be neglected. Therefore the 
temperature dependence of the SDT depends on the com- 
petition between the effect of the increase of the spin de- 
phasing due to the increase of the inhomogeneous broad- 
ening and the increase of the scattering with temperature 
(The latter provides additional spin dephasing channel) 
(Effect I) and the effect of the decrease of the spin de- 
phasing due to the counter effect from the increase of 
scattering which suppresses the inhomogeneous broaden- 
ing (Effect II). The results for impurity free case shown 



in Fig. 2(a) indicate that when T < 200 K and the total 
scattering is not so strong, Effect I is more important 
and hence the SDT decreases with T. When tempera- 
ture keeps increasing and the total scattering is further 
enhanced, the counter effect of the scattering to the inho- 
mogeneous broadening (Effect II) becomes more impor- 
tant and the SDT increases with T. Comparing with our 
previous works, 19 ' 21 one further finds that the absolute 
value of the SDT of LHs is one or two orders of magnitude 
smaller than the SDT of electrons. This can be easily un- 
derstood from the fact that the Rashba coefficients here 
are larger. 

Now we include the hole-impurity scattering with the 
impurity density = O.ITV^. As expected, when the to- 
tal scattering becomes stronger, the counter effect of the 
scattering to the inhomogeneous broadening takes the 
leading place easier and the SDT starts to increase with 
temperature earlier than the impurity-free case. When 
Ni = l.ONh, the total scattering is further enhanced. 
Now if one uses the hole- impurity scattering and the hole- 
phonon scattering to calculate the momentum relaxation 
time, and takes the lowest order of t" 1 ^) and tt LH (k) 
after expanding them over the function ^4j(#k) defined 
in Eq. (Bl), one gets the typical value of |J7 iH (k)|r p (k) 
at the average of k to be 0.11 at T = 100 K. It has al- 
ready entered the regime of strong scattering, and similar 
to the case of electrons when the linear part of the DP 
term is dominant, the SDT increases monotonically with 
temperature. 

When the well width becomes larger, the cubic term 
becomes more important. For example, when a = 7 nm, 

752 6i i^TOo7T§{ changes from —1.39 nm to —0.60 nm when 

the temperature changes from 100 K to 300 K. One can 
see from Table II that the cubic terms weighted by 
are dominant. The SDT in this case is plotted in Fig. 
2(a) as dashed curves for comparison. One can see that 
now the SDT decreases monotonically with temperature. 
It is because the increase of the inhomogeneous broaden- 
ing with temperature is much faster when the cubic term 
in the DP term is dominant. Therefore Effect I always 
surpasses Effect II with the increase of temperature and 
the SDT decreases monotonically with T. The same sit- 
uation happens in the case of HHs where there is only 
cubic term in the DP term. It is seen from Fig. 2(b) that 
the SDT decreases monotonically with temperature even 
when a = 5 nm. This is consistent to the electron case 
when the cubic term is dominant or the only term (bulk 
case) where the SDT also decreases monotonically with 
temperature. 17 ' 21 

One can find from the discussion above that the spin 
dephasing is a combined effect from the scattering and 
the inhomogeneous broadening due to the DP term. The 
inhomogeneous broadening induced spin dephasing 17 ' 18 
and the counter effect of the scattering to the inhomo- 
geneous broadening, are both very important and nei- 
ther can be neglected. Nevertheless, these effects arc 
either not or not fully accounted in the simplified model 
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which is based on the single-particle approximation Eq. 
(1). Furthermore, one should notice that the simplified 
model is based on the assumption of \fl\r p <C l, 6 which 
is not always satisfied for holes. To show the differences 
between the many-body approach and the earlier treat- 
ment, we also compare our results with those given by 
the simplified model which now reads: 
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in which 

r LH (k) 



,C dE k (f kXi+ - f k x,-)Tx(k) 
2 Jo°° dEk(fk\,+ — fk\,-) 



(18) 



+ T 3 , Lff ( 75 T) 2 ] , 



(19) 

r HH (k) = k°[n, HH ( 7 irr+Ts,HH(iIZ 7h ) 2 } ,(20) 



with 



/.A 



27T 



o- x (E k ,6)[l-cos{W)}d6 . (21) 



a\(E k ,0) stands for the scattering cross section of the 
hole-phonon and the hole-impurity scattering, and the 
expressions can be found in Eq. (B8). 

In the insets of Figs. 2(a) and (b), we plot the cor- 
responding SDTs of LHs and HHs from the simplified 
treatments in solid curves. SDTs from our many-body 
approach are plotted in dashed curves. One can see that 
both the curvatures and the absolute values are markedly 
different between the two treatments. The simplified 
treatment shows that the SDT of both HHs and LHs 
decreases monotonically with temperature regardless of 
the impurity densities. Moreover, when the impurity 
density increases, the SDT increases very fast. This is 
because the single-particle treatment totally ignores the 
fact that in the presence of the inhomogeneous broaden- 
ing, adding a new scattering means adding a new spin 
dephasing channel. It also does not treat the counter ef- 
fect of the scattering to the inhomogeneous broadening 
sufficiently. Moreover, it does not include the Coulomb 
scattering which we will show in the next subsection to 
be very important. By comparing the SDTs predicted by 
the two models, one can see that it is important to study 
the SDT of holes from the many-body approach. 



C. Effect of Coulomb scattering on SDT 

Now we turn to study the effect of the Coulomb scat- 
tering in p^\^a'\scatt to the SDT. It has been shown 
recently by Wu et al. that unlike the common belief 
that the Coulomb scattering cannot cause spin dephas- 
ing, in the presence of inhomogeneous broadening, it can 
also lead to spin dephasing 17 and for electrons in GaAs 
QWs, the Coulomb scattering is very important and can 
markedly increase the SDT. 19 ' 20 Glazov and Ivchenko 





FIG. 3: SDT vs. the scale coefficient of the DP term \ f° r 
(a) T = 120 K and (b) T = 300 K. The solid (dashed) curves 
are the results with (without) the Coulomb scattering. ■: 
Electrons; LHs; •: HHs. 



have also drawn the similar conclusion. 47 Since the spin- 
orbit coupling of hole system is much larger than that 
of electron one, it would be interesting to see how the 
Coulomb scattering can affect the SDT. Unlike the case 
of electron system, here we find that the hole-hole scatter- 
ing markedly reduces the SDT. This is consistent to the 
optical dephasing of semiconductors where the Coulomb 
scattering gives rise to a stronger optical dephasing. 17 ' 42 
In order to understand the difference between the hole 
system and the previous electron one, we plot in Fig. 3 
the SDT of the LH and the HH as function of a dimcn- 
sionless scale coefficient of the DP term \ at T = 120 
K and 300 K. Here x is introduced by hand in front 
of the DP term, i.e., x^ A (k) with \ = 1 correspond- 
ing to the case of the original DP term. The solid 
curves are for the case with both the hole-hole scatter- 
ing and the hole-phonon scattering and the dashed ones 
are for the case with the hole-phonon scattering only. 
It is pointed out here that notwithstanding the fact we 
sweep x through two orders of magnitude, experimen- 
tally the value of the Rashba coefficients can only be 
tuned within a small range by applying an external elec- 
tric field, 36,37 and can be determined by analyzing the 
Shubnikov-de Haas oscillations. 38 ' 39,40 ' 41 It is seen from 
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the figure that when x = 1, the SDTs of both the HH 
and the LH decrease when the Coulomb scattering is in- 
cluded. However, when one decreases the spin-orbit cou- 
pling by decreasing the scale coefficient x, one enters a 
regime where the Coulomb scattering increases the SDT. 
This is consistent with our previous observation that the 
competing effects of the inhomogeneous broadening and 
the scattering in different regimes. In the regime where 
the inhomogeneous broadening is weak (|$l A |r p <C 1), 
the hole-hole scattering mainly suppresses the inhomo- 
geneous broadening and consequently raises the SDT. 
In the regime where the inhomogeneous broadening is 
strong, (|J7 A |r p ^ 1), adding a new scattering provides 
a new spin dephasing channel and reduces the SDT. ft 
happens that for hole system, the spin-orbit coupling 36 
is within the strong inhomogeneous broadening regime. 
The same is true also for the optical dephasing where 
adding a new scattering also provides a new dephasing 
channel. 17 ' 42 One also finds from Fig. 3 that the SDT 
decreases with \. This can be easily understood because 
the spin dephasing becomes stronger when the spin-orbit 
coupling is larger. It is further noted that when the tem- 
perature rises, the effect from the Coulomb scattering is 
smaller. This is because the hole-phonon scattering be- 
comes more important with the increase of temperature. 

In order to compare the hole system with electron one, 
we add a scale coefficient \ in front of the electron Dres- 
selhaus term Eq. (2) and calculate the SDT of electrons 
in a QW (a = 5 nm) as function of X with and without 
the electron-electron Coulomb scattering. The results are 
plotted in the same figure for comparison. Similar to the 
case of holes, one finds that when the spin-orbit coupling 
is strong, the Coulomb scattering reduces the SDT also. 
It happens that the unsealed DP term (x = 1) is within 
the regime of weak inhomogeneous broadening. It is also 
seen from the figure that the SDT of holes is much smaller 
than that of electrons when X = 1- 



D. Impurity density dependence of SDT 

Now we turn to study how impurities affect the hole 
SDT. In Fig. 4 we plot the SDT of LHs as function of 
the impurity density Aj. It is seen from Fig. 4(a) that 
for LHs when X = 0.128 or X = 0.256, the SDT in- 
creases monotonically with the impurity density. Again 
the spin-orbit coupling here is in the regime where the in- 
homogeneous broadening is weak and the hole-impurity 
scattering mainly suppresses the inhomogeneous broad- 
ening and raises the SDT. When the scale coefficient 
X = 0.512 and Ni = 0, \fl LH (k)\T p (k) is close to 1 at 
the average of k: [\fl LH (k)\T p (k)]\ k={k) = 0.32. The ef- 
fect that adding a new scattering provides a new spin 
dephasing channel becomes dominant and the SDT first 
decreases with the impurity density. This is similar to 
the effect of the Coulomb scattering discussed above 
while the spin-orbit coupling is large. However, with 
the increase of the impurity density, r p gets smaller and 




FIG. 4: SDT vs. the impurity density Ni/N h with different 
scale coefficients of the DP term \ at T — 120 K for (a) LHs 
and (b) HHs. ■: x = 0.128; A: X = 0.256; •: X = 0.512; + 
X = 1-024. 



\n LH (k)\T p (k) gets smaller again. When N { = 0.1N h , 
[\n LH {k)\T p (k)]\ k={k) = 0.16, the SDT reaches a mini- 
mum. Further increasing the impurity density, one enters 
the strong scattering regime ([\fl LH (k)\T p (k)]\ k= ^ = 
0.049 « 1 at JV t = l.0N h ), and the SDT keeps increas- 
ing with N in this regime. When the scale coefficient 
X = 1.024, | O (k) | becomes even larger and the minimum 
of the SDT occurs at larger ATj. The similar is also true 
for the case of HHs in Fig. 4(b). It is noted that in re- 
ality x is around 1 and the SDT will first decrease then 
increase with the impurity density. This is totally dif- 
ferent from the electron case and also different from the 
prediction of the single-particle approach where the SDT 
always increases with the impurity density. 

We further investigate the impurity density depen- 
dence of the SDT at different temperatures. Here the 
scale coefficient X is fixed to be 1. In Fig. 5 one finds 
that when the temperature is low, the SDT first decreases 
with Ni as then the total scattering is weak, and then in- 
creases with it after the SDT reaches a minimum. When 
the temperature gets higher, the SDT always increases 
with A^j as then |fi A |r p <C 1 can always be satisfied and 
the hole-impurity scattering mainly suppresses the inho- 
mogeneous broadening. 
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FIG. 5: SDT vs. the impurity density at different tempera- 
tures for (a) LHs and (b) HHs. \ = 1 



E. Hole density dependence of SDT 

Finally we investigate the hole density dependence of 
the SDT at different temperatures and well widthcs. 
Here the hole- impurity scattering is excluded and x = 1 . 
In Fig. 6(a) we plot the SDT of LHs as function of the LH 
density with a = 5 nm. The SDT decreases with the hole 
density when the temperature is low but increases with 
it when the temperature is high enough. To understand 
this result, we first analyze the Rashba term [Eqs. (4-9)]. 

When T = 300 K, one finds j®!f l E z mojT?l changes from 
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fcf) 

—2.58 nm to —1.46 nm when Nh changes from 5 x 10 
cm~ 2 to 4 x 10 11 cm -2 , and the absolute value becomes 
even larger when T = 100 K. Therefore, it can be seen 
from Table II that both jf l 3 61 and 7 ff 1 are smaller than 

— 752 6 '|||y, and the linear terms in Eqs. (7) and (8) are 

dominant. Moreover, again |f2 Lff |r p is slightly smaller 
than 1. Similar to the case in Sec. Ill B when the linear 
Rashba term is dominant, the hole density Nh influences 
the spin dephasing through two competing effects: Effect 
I: The increase of the spin dephasing due to the increase 
of the inhomogeneous broadening with Nh as holes are 
populated at higher fc-states at high hole density; and 
due to the increase of the scattering which provides ad- 



100 K 
140 K 
180 K 

a = 5 nm 220 K 

300 K 




N h (10 11 cm" 2 ) 



FIG. 6: SDT vs. the hole density at different temperatures, 
(a): LHs with a = 5 nm; (b): LHs with a = 7 nm; (c): HHs 
with a — 5 nm. 



ditional spin dephasing channel. Effect II: The decrease 
of the spin dephasing due to the counter effect of the in- 
creased scattering which suppresses the inhomogeneous 
broadening. The results shown in Fig. 5(a) indicate that 
when T < 220 K and the total scattering is not so strong, 
Effect I is more important and the SDT decreases with 
Nh. When the hole density keeps increasing and the to- 
tal scattering is further enhanced, Effect II becomes more 
important and the SDT increases with N h . 

We further plot the SDT of LHs with a = 7 nm in Fig. 
6(b) where the cubic terms become more important. For 
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example, when T = 300 K, 7?n 6i E 'zTOottIt changes from 

— 1.06 nm to —0.60 nm when Nh changes from 5 x 10 10 
cm~ 2 to 4 x 10 11 cm -2 . One can see from Table II that 
the cubic terms weighted by 753 s ' are dominant. Similar 
to the case in Sec. Ill B, when the cubic Rashba term is 
dominant, the increase of the inhomogeneous broadening 
with hole density is much faster than the counter effect of 
the scattering and consequently Effect I always surpasses 
Effect II with the increase of hole density. As expected, 
one finds that the SDT decreases monotonically with Nh- 
The same is true for HHs in Fig. 6(c) where the Rashba 
term [Eqs. (4-6)] includes only the cubic one. 



IV. CONCLUSION 

In conclusion, we have performed a systematic micro- 
scopic many-body investigation on the hole spin dephas- 
ing of p-type GaAs QWs of small well width where the 
HH and LH bands are well separated, by constructing a 
set of kinetic spin Bloch equations based on the noncqui- 
librium Green function method. We included the mag- 
netic field, the Rashba spin-orbit coupling and all spin 
conserving scattering such as the hole-phonon, the hole- 
nonmagnetic impurity and the hole-hole scattering. By 
numerically solving the kinetic equations, we obtained 
the time evolution of the distribution functions and the 
spin coherence of holes. The SDT is calculated from the 
slope of the envelope of the incoherently summed spin co- 
herence. Differing from earlier studies on spin dephasing 
based on the single-particle approach which only includes 
the lowest-order elastic scattering and the anisotropy 
from — |J~2(k)| and +|f2(k)|, this approach takes full ac- 
count of the inhomogeneous broadening from different 
k-states of the Rashba term as well as the effect of all 
the scattering. Furthermore, this approach is valid re- 
gardless of the strength of scattering whereas the earlier 
single-particle approach is valid only when the scattering 
is strong enough, i.e., |J7|r p -C 1. Using this many-body 
approach, we studied in detail how the hole spin dephas- 
ing are affected by temperature, the hole-hole Coulomb 
scattering, the impurity and the hole densities. 

We showed that the spin dephasing is mainly affected 
by two effects: The inhomogeneous broadening and the 
scattering. Any effect that increases the inhomogeneous 
broadening tends to reduce the SDT. However, the ef- 
fect of scattering on the spin dephasing is different when 
\n\r p < 1 and \n\r p 1: When |fi|r p < 1 and there- 
fore the scattering is strong in comparison to the DP 
term, the counter effect of the scattering to the inho- 
mogeneous broadening is important. In this regime, the 
scattering tends to drive carriers to a more homogeneous 
state in k-space and consequently reduces the inhomoge- 
neous broadening. This tends to increase the SDT. When 
\n\r p gl 1, the scattering is weak in comparison to the DP 
term (inhomogeneous broadening) and the counter effect 
can be neglected, adding a new scattering provides an ad- 



ditional dephasing channel. In this regime, the counter 
effect of the scattering to the inhomogeneous broadening 
can be ignored and the scattering reduces the SDT. All 
the factors, such as temperature, well width, impurity 
density and hole density, can affect the inhomogeneous 
broadening and the scattering and therefore influence the 
SDT. 

The temperature affects the SDT in two ways: On the 
one hand, the increase of the temperature drives holes to 
higher k-states, and leads to a stronger inhomogeneous 
broadening. On the other hand, the scattering is en- 
hanced with the increase of the temperature. When the 
linear Rashba term is dominant, such as LHs with a = 5 
nm at the hole density in our investigation, it is shown 
that the SDT decreases with T when the temperature is 
low and the impurity density is small. This can be un- 
derstood as it is in the regime where | | and l/r p are 
comparable and the increase of the spin dephasing due 
to the increase of the inhomogeneous broadening and the 
increase of the spin dephasing channel by the increase of 
the scattering with temperature are dominant. When the 
temperature keeps increasing so that the scattering be- 
comes stronger or when the impurity density is high, the 
SDT increases with T when the system enters the regime 
where \CI lh \t p <C 1 and the counter effect of the scat- 
tering to the inhomogeneous broadening becomes dom- 
inant. When the cubic Rashba term is dominant (such 
as LHs with a = 7 nm in our investigation) or is the 
only term (such as HHs), the SDT decreases monotoni- 
cally with temperature as the increase of the inhomoge- 
neous broadening with temperature is much faster than 
the increase of scattering. These results are quite differ- 
ent from the case of electrons where the spin-orbit cou- 
pling is within the regime of weak inhomogeneous broad- 
ening (|J~i|T p <C 1) and the SDT increases monotonically 
with temperature when the linear DP term is dominant. 
We also compared the SDTs predicated by our many- 
body approach with the results of the earlier simplified 
treatment, and showed that the simplified treatment is 
inadequate in studying the hole spin dephasing. 

The hole density also influences the inhomogeneous 
broadening and the scattering simultaneously. Similar 
to the case of temperature dependence, it is shown that 
for the LHs with a = 5 nm where the linear Rashba 
term is dominant, the SDT decreases with Nh when the 
temperature is low because it is in the regime of strong 
inhomogeneous broadening, and increases with Nh when 
the temperature is higher and the inhomogeneous broad- 
ening is weak. For LHs with a = 7 nm or HHs where the 
cubic Rashba term is the leading/only term, the SDT 
decreases monotonically with Nh- 

We further showed that the Coulomb scattering con- 
tributes markedly to the SDT. When the inhomogeneous 
broadening is stronger than the scattering, the Coulomb 
scattering enhances the spin dephasing. Otherwise, it re- 
duces the spin dephasing. In the earlier single-particle 
treatment, the Coulomb scattering was considered to be 
unable to cause spin dephasing. 
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In the calculation, the magnetic field in the Voigt con- 
figuration is taken to be 4 T. We found that for hole sys- 
tem, the magnetic field dependence is marginal as the 
Rashba term is very large. In this investigation, the 
Elliott- Yafet mechanism 48 is not included. A full mi- 
croscopic many-body treatment of the this mechanism 
is much more complicated than the DP mechanism and 
will be published elsewhere. Up till now there is no ex- 
perimental investigation on the SDT for holes in (001) 
QWs. Experiments such as spin-echo measurements 49 
and time-resolved Faraday rotation measurements 50 can 
be used to measure the SDT. 
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APPENDIX A: EFFECT OF SCATTERING TO 
SDT 

In Sec. Ill B we pointed out that when the scattering is 
weak in comparison to the DP effective field (inhomogc- 
neous broadening) , the main effect of scattering is to add 
an additional dephasing channel. Therefore the scatter- 
ing here reduces the SDT. To demonstrate this effect, we 
now study the spin dephasing in the limiting case of no 
scattering included in the calculation. In Fig. 7 we show 
the temporal evolution of the incoherently summed spin 
coherence plh = X^k \PkLH{t)\ for LHs with total LH 
density N h = 4 x 10 11 cm" 2 and T = 300 K. The coher- 
ently summed spin coherence p' LH — |^ k pwLn(t)\ is also 
plotted for comparison. One can see that the amplitude 
of p' LH oscillates and decays to zero very quickly due to 
the interference caused by the momentum dependence of 
the DP term. However, the incoherently summed spin 
coherence plh does not decay, which means that there 
is no irreversible dephasing. 16,42,44 This is consistent to 
the fact of no scattering as there is no dissipation pro- 
cess here. By adding a scattering, one introduces the 
dissipation into the system which causes an irreversible 
dephasing. This can be seen in the same figure where the 
incoherently summed spin coherence plh is plotted with 
the impurity density Ni — O.OIA^. One finds that the 
Plh now decays with time, although much slower than 
the one in Fig. 1 where all the scattering is included. 
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FIG. 7: Temporal evolution of the incoherently summed spin 
coherence pLu{t) (solid curve) and the coherently summed 
spin coherence Plh^) (Dotted curve) without any scatter- 
ing and the temporal evolution of the incoherently summed 
spin coherence pLH{t) with only hole-impurity scattering with 
Ni = 0MN h (Chained Curve). 



APPENDIX B: A SIMPLIFIED ANALYTICAL 
ANALYSIS OF SDT 

In Sec. Ill B we pointed out that the ratio of the DP 
term to the scattering rate determines the way how the 
scattering affects the spin dephasing. To reveal part of 
this effect analytically, we now study a much simplified 
case with only the hole-impurity scattering and the dom- 
inant part of the DP term, i.e., for LHs we include only 
the linear Rashba term weighted by and for HHs 
only the cubic Rashba term weighted by 754 7h . Further- 
more, we will also neglect the inhomogeneous broadening 
later. 

First we consider the HH case. We expand 2x2 density 
matrix p^nn as follows: 



PkHH = ^PHHd(k) A l(^) 



(Bl) 



with Ai(6k) = -±=e iWk . The coherent terms of the 
kinetic spin Bloch equations [Eqs. (14) and (15)] can 
be written into the matrix commutator as pkHH\coh = 
t [H^ H (k), PKHH ] with H» H (k) = \a- n" H (k). It 
is noted that here we neglected the Coulomb Hartree- 
Fock term. Furthermore, we expand H^ H as H^ H (k) = 
J2i H^ t H (k)Ai(6k). Therefore the coherent term in the 
matrix form can be written as: 



P\t,HH\coh 



E 

Mi 



[H^ h (k), PHHM (k)]Me k ). 



(B2) 

With only the dominant part of the DP term (term 
weighted by 754 7/1 ) included, H^ H (k) is expanded as: 



ttHH 

H s-3 - 

ttHH 



-iS-y 54 
, 



(B3) 
(B4) 
(B5) 
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in which S = \{p x — ia v ). Substituting Eqs. (B3-B5) 
into Eq. (B2), one obtains 

p k HH\coh - -^=^ 7h E z klJ2([S, PHH , l+3 (k)} 

-[S^, P HH,i-3(k)])Ai(9k) • (B6) 
Similarly one can expand the scattering term as: 

pkHHlacat = Y, PHH,l{k)Uf{k)A l {6) , (B7) 



in which 



for the elastic scattering with 



0f(»).arw ( fjf. 



d9 

(27T) 



^(l-cosW). (B8) 



Here q{9) = ^J2k 2 {\ - cos 9), and U 2 (e) = 

£„{47re 2 /M<Z 2 (0) + <£)]}W&)| 2 is the hole- 
impurity scattering matrix element and \I(iq z )\ 2 = 
7T 2 sin 2 y/[y 2 (y 2 - 7r 2 ) 2 ] with y = q z a/2 is the form 
factor. 

Now we can expand the spin Bloch equations Eq. (13) 
into pHH,i{k) as follows: 

PHH.l(k) + ^f h E z k\^= ([S, PHH.l+m\ 

- [S\pHH.i- 3 (k)}) = -Uf{k) PH H,i{k) .(B9) 

In order to find the solution, we multiply er = (a x , cr y , cr z ) 
to both sides of this equation and calculate the trace. By 
defining the HH "spin" vector to be Tr (pHH.i(k)cr) = 
&HH,i{k), one can rewrite Eq. (B9) into 

S HH ,i(k) +lli 7h E z kl-^=(T ■ S HH , l+3 (k) 

V Z7T v 



- P ■ S HH ,i-3(k)) = -C/ ; 2 (fc)S Hff , ; (fc) ,(B10) 

thanks to the relation Tr([5, pHH,i-i{k)]cr) = 
Tr(p H H,l-i(k)[(T,S\). In Eq. (BIO) the tensor T 
reads 



T = 



I 1 \ 

-i 
[-1 i J 



(Bll) 



One finds that SnH,i(k) is only related to S##,;±3(fc). 
By considering the lowest orders with I = 0, and ±3 and 
defining Shh {k) = (S HH - 3 (k), S H H,o( k )> s ra,3W) T , 
Eq. (BIO) can be written as: 



(k) + ( 



7h7h EM 



754 



2tt 



-Q+U)S HH {k) = , (B12) 



W = 



u!(k) 




(B13) 



(B14) 



In Eq. (B14) we have used the relations U 2 (k) = U 2 3 (k) 
and [7 2 (fc) =0. 

Now one can solve Eq. (B12) analytically. To reveal the 
main characteristic analytically, we make the assumption 
that the spin relaxation/dephasing occurs mainly around 
the Fermi surface and 6\s_ f = 0. By doing so, one throws 
away the interference between different k-states [except 
for the states with 0k f (— here), and 9-^ } ± 2n/3], and 
therefore the inhomogeneous broadening. Then the HH 
spin Shh = E { s hh i > s hh n S HHi) has the form: 

/=-3,0,3 

S HH = , (B15) 

o -iti(x+v / 5 5zr T6) 



y/x 2 - 16 
e -±ti(a:+Vz 2 -16) 



x e 



tiVa; 2 -16 



1) 



2Vs 2 - 16 

+ Vx 2 - 16(e tlVx '- 16 + i)]S , (B17) 



u- u Ui(kf) 

m wmcn a; — 



/2tt 



^7h7hp 1.3 
'54 ^z^f 



is proportional to the ratio of the 



scattering rate to the DP term, 



t and 5n 



represents the initial spin polarization along the z-axis. 
One can see from Eq. (B16) that when x < 4, the SDT is 
proportional to 1/x and decreases with x whereas when 
x > 4, the SDT is proportional to l/{x— \J x 2 — 16) and 
increases with x. This result indicates that the scattering 
reduces the SDT when the scattering is weak in compari- 
son to the DP effective field but increases the SDT when 
the scattering is strong in comparison to the DP effective 
field. Moreover, when U 3 = 0, i.e., there is no scatter- 
ing, x — and consequently there is no spin dephasing. 
This is consistent with the numerical result presented in 
Appendix A. 

Similarly one can derive the equation for spin of LHs 
Slh with only the linear part of the Rashba term in- 
cluded. One gets the same equations Eqs. (B15-B17) but 

with x - ltrEz[kl)kf and H - ^- t. 

These results coincide qualitatively with the results 
shown in Fig. 4: The SDT first decreases then increases 
with the hole-impurity scattering when the spin-orbit 
coupling is strong; but increases monotonically with the 
scattering when the spin-orbit coupling is weak. Further- 
more, by making the approximation that the hole-phonon 
scattering is also an elastic scattering and by including 
the hole-phonon scattering in the scattering term, U 2 in 
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Eq. (B8) can be modified as: 

Uf{k) = 2^fj JtL [Nl UZ + (l + 2N q ) 9 i} 

x (1 -cosW) . (B18) 

Here 3g 2 -E{2 7 r e 2 fi L o/[(g 2 + g z 2 )]}(l/ K co-l/ K o)|/fe)| 2 

is the holc-phonon interaction matrix element, and N q — 
l/[exp(il LO /fc B T) — 1] is the Bose distribution of the LO 



phonon. Then the results also coincide qualitatively with 
those in Fig. 5(b): the SDT first decreases then increases 
with the hole-impurity scattering when the total scatter- 
ing is weak, but always increases with scattering when 
the total scattering is strong. Finally we point out that 
as we do not include the inhomogeneous broadening and 
all the scattering in this simplified model, many other 
features predicted in the text cannot be obtained here. 
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